Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 27 December 2012 (MN WgX style file v2.2) 



A new multi-field determination of the galaxy luminosity function at 
z = 7 — 9 incorporating the 2012 Hubble Ultra Deep Field imaging 

R. J. McLure 1 *, J. S. Dunlop 1 , R.A.A. Bowler 1 , E. Curtis-Lake 1 , M. Schenker 2 , 
R.S. Ellis 2 , B.E. Robertson 3 , A. M. Koekemoer 4 , A.B. Rogers 1 , Y. Ono 5 , M. Ouchi 5 , 
S. Chariot 6 , V. Wild 7 , D.R Stark 3 , S.R. Furlanetto 8 , M. Cirasuolo 1 , T.A. Targett 1 

1 SUPA'\, Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh EH9 3HJ 

2 Department of Astrophysics, California Institute of Technology, MS 249-17, Pasadena, CA 91125, USA 

3 Department of Astronomy and Steward Observatory, University of Arizona, Tucson AZ 85721, USA 
4 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 

° Institute for Cosmic Ray Research, University of Tokyo, Kashiwa City, CHiba 277-8582, Japan 
e UPMC-CNRS, UMR7095, Institut d'Astrophysique de Paris, F-75014, Paris, France 

7 School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews, KY16 9SS 

8 Department of Physics & Astronomy, University of California, Los Angeles, CA 90095, USA 



27 December 2012 



ABSTRACT 

We present a new determination of the UV galaxy luminosity function (LF) at redshift z ~ 7 
and z ~ 8, and a first estimate at z ~ 9. An accurate determination of the form and evo- 
lution of the galaxy LF during this era is of key importance for improving our knowledge 
of the earliest phases of galaxy evolution and the process of cosmic reionization. Our analy- 
sis exploits to the full the new, deepest WFC3/IR imaging from our Hubble Space Telescope 
(HST) Ultra Deep Field 2012 (UDF12) campaign, with dynamic range provided by includ- 
ing a new and consistent analysis of all appropriate, shallower/wider-area HST survey data. 
Our new measurement of the evolving LF at z ~ 7 — 8 is based on a final catalogue of 
~ 600 galaxies, and involves a step-wise maximum likelihood determination based on the 
photometric -redshift probability distribution for each object; this approach makes full use of 
the 11-band imaging now available in the Hubble Ultra Deep Field (HUDF), including the 
new UDF12 F140W data, and the latest Spitzer IRAC imaging. The final result is a deter- 
mination of the z ~ 7 LF extending down to UV absolute magnitudes M1500 = —16.75 
(AB mag), and the z ~ 8 LF down to A/1500 = —17.00. Fitting a Schechter function, we 
find M* 5Q0 = — 19.90±g;H, lo S^* = -2-96^;^, and a faint-end slope a = -1.90±8;^ 
at z ~ 7, and M* 5m = -20.12t°;g, log0* = — 3.35±g;||, a = — 2.02±g;|| at z ~ 8. 
These results strengthen previous suggestions that the evolution at z > 7 appears more akin 
to 'density evolution' than the apparent 'luminosity evolution' seen at z ~ 5 — 7. We also pro- 
vide the first meaningful information on the LF at z ~ 9, explore alternative extrapolations 
to higher redshifts, and consider the implications for the early evolution of UV luminosity 
density. Finally, we provide catalogues (including derived z p h a u A/1500 an d photometry) for 
the 100 most robust z ~ 6.5 — 11.9 galaxies in the HUDF used in this analysis. We briefly 
discuss our results in the context of earlier work and the results derived from an independent 
analysis of the UDF12 data based on colour-colour selection (Schenker et al. 2013). 

Key words: galaxies: high-redshift - galaxies: evolution - galaxies: formation 



1 INTRODUCTION 

The advent of deep near-infrared imaging, in particular with Wide 
Field Camera 3 (WFC3/IR) on the Hubble Space Telescope (HST) 
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has now enabled the discovery and study of galaxies to be ex- 
tended to redshifts z ~ 6.5 — 10, into the first billion years of 
cosmic history (see Dunlop 2012 for a review). This work is of 
fundamental importance for improving our understanding of the 
formation and growth of the early generations of galaxies, and test- 
ing the predictions of the latest galaxy-formation simulations. It is 
also of interest for establishing whether these galaxies reionized 
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the Universe (e.g. Robertson et al. 2010; Finkelstein et al. 2012; 
Kuhlen & Faucher-Giguere 2012), and, if so, providing more de- 
tailed information on how reionization proceeded (as compared to 
the integrated constraints on 'instantaneous' reionization provided 
by current measurements of microwave background polarization - 
z re i 0n ~ 10.6 ± 1.2; Komatsu et al. 201 1). 

Because galaxies at z ~ 7 are so faint, it is hard to gain de- 
tailed physical information on the properties of individual objects, 
and indeed only a handful of spectroscopic redshifts have been es- 
tablished on the basis of Lyman-a emission at z ~ 7 (the current 
record holder is at 2 = 7.213; Ono et al. 2012). Attention has 
thus (sensibly) focussed on population statistics, helped by the fact 
that significant samples of photometrically-selected galaxies can 
now be assembled at these redshifts due to the presence of a strong 
Lyman-break at A res t ~ 1216 A in their spectral energy distribu- 
tions (SEDs). This is caused by near-complete absorption by neu- 
tral hydrogen gas along the line-of-sight at z > 6.5 (Fan et al. 2006; 
Mortlock et al. 2011), making 'Lyman-break' galaxy (LBG) selec- 
tion in principle straightforward at these redshifts, given adequate 
data. 

The first and most important population measurement which 
is usually attempted once a significant sample of galaxies is avail- 
able at a given redshift is a determination of the luminosity function 
(LF); i.e. the comoving number density of galaxies as a function 
of luminosity (= absolute magnitude). Prior to the 2009 installa- 
tion of WFC3/IR the availability of deep zsso and 4775 imaging 
from the Advanced Camera for Surveys (ACS) on HST enabled the 
UV (Xrest — 1500 A) LF to be established for faint galaxies out 
to 2 ~ 6 by Bouwens et al. (2007). In a complementary effort 
based on the new availability of degree-scale red/infrared imag- 
ing from ground-based telescopes, the bright end of the UV LF 
was also measured out to z ~ 6 by McLure et al. (2009), who 
demonstrated that a combined analysis yielded a consistent result. 
Both of these studies exploited the available photometry to estab- 
lish the presence/location of the afore-mentioned Lyman-break, but 
whereas Bouwens et al. (2007) continued with the simple and well- 
established two-colour selection technique, McLure et al. (2009) 
used SED fitting with evolutonary synthesis models (e.g. Bruzual 
& Chariot 2003) to derive photometric redshifts based on all of the 
available multi-band photometry. 

These two alternative approaches to galaxy selection have now 
both been exploited to explore the form of the LF at higher red- 
shifts. Specifically, with the new first-epoch deep WFC3/IR data 
provided by the UDF09 program (GO 11563, PI: Illingworth), 
Oesch et al. (2010a) and McLure et al. (2010) produced alterna- 
tive (but again consistent) determinations of the z ~ 7 galaxy UV 
LF. The SED fitting approach was also exploited by Finkelstein et 
al. (2010), while colour-colour selection has since been re-applied 
by Bouwens et al. (201 la) at 2 ~ 7 and z ~ 8, to the final UDF09 
dataset. Colour-colour selection has also recently been applied in 
attempts to constrain the brighter end of the LF at z ~ 8 by Bradley 
et al. (2012) and Oesch et al. (2012), to the BoRG^and CANDELS 
(Grogin et al. 201 lQdatsets respectively. Meanwhile, SED fitting 
has been applied to the UDF09 and CANDELS GOODS-South 
data by Finkelstein et al. (2012) in order to derive a new estimate 
of the evolving UV luminosity density, and by Bowler et al. (2012) 



in the search for brighter z ~ 7 galaxies in the early Ultra VISTA 
data (McCracken et al. 20120 

It is important to note that, given only three-filter data (i.e. 
two-wavebands above a putative Lyman-break and one below) 
colour-colour selection and SED fitting are essentially equivalent. 
However, as the number of useful wavebands expands, it is clear 
that SED fitting makes more complete and consistent use of the 
available data. This has, to some extent, been recognized by the 
adoption of additional somewhat ad-hoc criteria to colour-colour 
selection, in an attempt to 'factor-in' the extra information provided 
through other filters (e.g. the rejection of objects which which show 
more that one > 1.5-cr detection in bluer bands and the compu- 
tation of a separate Xopticai by Bouwens et al. 2011a). However, 
SED fitting clearly deals with all detections and non-detections in a 
more straightforward and consistent manner, and has the additional 
benefit of providing actual redshift estimates with confidence in- 
tervals (and indeed can provide a redshift probability distribution 
for each object, albeit this depends somewhat on adopted priors; 
McLure et al. 201 1). Finally, SED fitting also more clearly exposes 
the nature of potential interlopers (such as dusty red galaxies, post- 
starburst objects with strong Balmer breaks, and dwarf stars in our 
own galaxy) and provides clearer information on what data needs 
to be improved to eliminate them (e.g. Dunlop 2012). Neverthe- 
less, whatever the selection technique, careful simulation work is 
required to quantify selection bias, completeness and contamina- 
tion in any determination of the evolving galaxy LF in the young 
universe. 

In an attempt to make further progress, and in particular to 
extend the study of galaxies both to higher redshifts (2 > 8) and 
lower luminosities (at 2 ~ 7 — 8) we have recently completed a 
new programme of even deeper near-infrared imaging in the Hub- 
ble Ultra Deep Field (HUDF; Beckwith et al. 2006) with WFC3/IR 
on HST. This new imaging campaign was completed in Septem- 
ber 2012 (GO 12498, PI: Ellis, hereafter UDF12) and when com- 
bined with the existing UDF09 data provides the deepest ever near- 
infrared images of the sky. The data have now been reduced and 
released to the public (Koekemoer et al. 2013) through the team 
websit^] Key elements of our observing strategy were the delivery 
of extremely deep Y105 data to more robustly identify galaxies at 
2 ~ 8 and higher redshifts, and the addition of imaging through 
a new filter previously unexploited in the HUDF, J140, both to en- 
able galaxy discovery to be pushed beyond 2 ~ 8.5, and to en- 
able more accurate SEDs to be determined for galaxies at 2 ~ 7 
and z ~ 8. The final UDF12+UDF09 combined dataset reaches 
the planned 5-cr detection limits of Y105 = 30.0, J125 = 29.5, 
J140 = 29.5, ifieo = 29.5 (in apertures of diameter 0.40", 0.44", 
0.47", 0.50" respectively, sampling 70% of point-source flux den- 
sity in each waveband). The results of our search for galaxies at 
2 > 8.5 have already been reported by Ellis et al. (2013), while the 
new deep multi-band data have now also been exploited by Dun- 
lop et al. (2013) in a new determination of the UV spectral slopes 
of galaxies at 2 ~ 7 — 9 (with consequent implications for their 
stellar populations). Most recently, a new determination of galaxy 
sizes based on the UDF12 dataset has been completed by Ono et al. 
(2013). 

In this paper we focus on utilising the UDF12 dataset, along 
with the ever-growing shallower WFC3/IR imaging over wider ar- 
eas, to undertake a new determination of the galaxy UV LF at 2 ~ 7 
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and z ~ 8. Crucially the new ultra-deep imaging in the HUDF im- 
proves our ability to probe the faint end of the LF, better sampling 
the population of numerous faint galaxies (h'hsoo > —18) which 
likely dominate the UV luminosity density and hence drive reion- 
ization. A key goal, therefore, is to better establish the faint-end 
slope, q, on which extrapolations to even fainter (as yet unobserv- 
able) luminosities have to be based. However, simply increasing 
the depth of the deepest field does not yield significantly better es- 
timates of a unless the degeneracies between the Schechter func- 
tion parameters (M* , <fi* , and a) can be minimized (e.g. Bouwens 
et al. 2011a; Dunlop 2012). This requires maximising the usable 
dynamic range in UV luminosity, to properly constrain the shape 
of the LF. Thus, to best determine the z ~ 7 and z ~ 8 LF, we 
have analysed the new HUDF 12 data in combination with the pro- 
gressively shallower WFC3/IR survey data provided by the UDF09 
parallel fields, the Early Release Science (ERS) data in GOODS- 
South, the CANDELS data in the remainder of GOODS-South 
and the UDS field, and all of the parallel BoRG data obtained by 
Sept 2012. In each of these fields, (including the HUDF) we have 
also utilised the associated HST ACS imaging, Spitzer IRAC data 
(deconfused with the WFC3/IR Hi6o imaging), an d ground-based 
near-infrared/optical data where appropriate. 

Given that the selection functions and associated simulations 
are different, and to facilitate comparison with other work, our team 
has also undertaken a parallel, and completely independent deter- 
mination of the LF at z ~ 7 and z ~ 8 based on 'traditional' drop- 
out colour-colour selection. The results from this are presented in 
Schenker et al. (2013) but are also summarised in this present pa- 
per for ease of comparison with the SED-fitting technique results 
derived here. 

The remainder of this paper is structured as follows. In Sec- 
tion 2 we give full details of the datasets utilised in this new study, 
and explain how we selected galaxy catalogues before refining the 
samples to contain only plausible high-redshift galaxy candidates. 
We also describe the simulations undertaken to establish complete- 
ness and contamination corrections in each of the individual survey 
fields, simulations which are crucial for a robust determination of 
the LF from such a complex multi-field dataset. Next, in Section 3, 
we describe how we chose to determine the LF, adopting as our pri- 
mary technique the non-parametric step-wise maximum likelihood 
method (SWML), but also applying parametric maximum likeli- 
hood fitting to explore Schechter-function representations of the 
LF. We then present the results of our analysis in Section 4, pro- 
viding our best measurements of the LF at z ~ 7, z ~ 8 and 
z ~ 9, and briefly exploring the implied evolution of the LF with 
redshift. Here we also compare our results with the independent 
UDF12 analysis of Schenker et al. (2013), and discuss our derived 
LF parameters (with associated improved confidence intervals) in 
the context of the the results deduced by Bouwens et al. (2011a), 
Bradley et al. (2012) and Oesch et al. (2012) prior to UDF12. In 
Section 5 we proceed to explore the implications of our results 
for the evolution of the LF out to even higher redshifts, and de- 
rive the implied evolution of UV luminosity density as a function 
of redshift (a key measurement for tracking the likely progess of 
reionization). Finally, we present a summary of our conclusions 
in Section 6. Throughout the paper we will refer to the follow- 
ing HST ACS+WFC3/IR filters: F435W, F600LP, F606W, F775W, 
F814W, F850LP, F098M, F105W, F125W, F140W & F160W 

as -B435, V(300, V606, *775, 1814, 2850, Yogs, Yl05, Jl25, Jl40 & -Hl60 

respectively. All magnitudes are quoted in the AB system (Oke 
1974; Oke & Gunn 1983) and all cosmological calculations assume 
fi Q = 0.3, = 0.7 and H = 70 kms _1 Mpc _1 . 



2 DATA 

In this section we provide a summary of the basic properties of the 
datasets which we have utilised in this study to measure the high- 
redshift galaxy luminosity function. In addition, we also provide 
details of the methods adopted to derive accurate image depth infor- 
mation, object photometry and reliable catalogues of high-redshift 
galaxy candidates. 

2.1 Survey fields 

The datasets analysed in this paper form a 'wedding-cake' struc- 
ture ranging from the ultra-deep UDF12 observations covering an 
area of only ~ 4.5 arcmin 2 to wider-area WFC3/IR survey data 
covering several hundred arcmin 2 . Below we provide the basic ob- 
servational details of each dataset in turn. 

2.1.1 The UDF12 

The dataset which plays the pivotal role in constraining the faint- 
end of the galaxy luminosity function at z > 7 and provides the 
primary motivation for this paper is the new UDF12 WFC3/IR 
multi-band imaging of the Hubble Ultra Deep Field (GO- 12498, 
RI. Ellis). The UDF12 observing campaign acquired 128 orbits of 
WFC3/IR integration time targeting the HUDF, all of which were 
obtained between 4th August and 16th September 20122 As dis- 
cussed in the introduction, the primary motivation for the UDF12 
observing campaign was to improve our knowledge of number den- 
sities and spectral properties of the ultra-faint galaxy population at 
z = 7 — 8 and to provide the first robust census of the z > 8.5 
galaxy population. 

In order to achieve these aims the bulk of the UDF12 orbits 
were invested in quadrupling the HUDF integration time in the cru- 
cial Y105 filter and providing ultra-deep imaging in the J140 filter, 
which had not been employed in previous HUDF imaging cam- 
paigns. The 128 orbits awarded to UDF12 were allocated as fol- 
lows: 72 orbits in the Y105 filter, 30 orbits in the J140 filter and 26 
orbits in H160. In combination with the data provided by the previ- 
ous UDF09 observing campaign (GO-11563, RI. Illingworth), the 
total orbit allocation of dedicated WFC3/IR data in the HUDF now 
stands at: 96 orbits in Y105, 34 orbits in J125, 30 orbits in J140 and 
79 orbits in H±eo- The depths of the available data in the HUDF 
(and the other survey fields analysed in this study) are provided in 
Table 1. 

2.1.2 The HUDF09 parallel fields 

To increase our ability to constrain the faint-end of the high-redshift 
luminosity function we have also utilised the WFC3/IR imaging 
available in the two parallel fields of the HUDF09 imaging cam- 
paign, hereafter HUDF09-1 and HUDF09-2. Although substan- 
tially shallower than the data available in the HUDF itself, the 
HUDF09 parallel fields consist of 33 and 48 WFC3/IR orbits re- 
spectively (spread between the Y105, J 125 & ^ieo filters) and pro- 
vide crucial leverage for constraining the z > 7 galaxy luminos- 
ity function ~ 1.5 magnitudes brighter than the ultimate depth 
achieved by UDF12. For the purposes of this analysis we have 
utilised our own reduction of the HUDF09 WFC3/IR imaging in 

5 The entire reduced UDF12 dataset is publicly available on the following 
web-site : http ://archive . stsci . edu/prepds/hudf 1 21 
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Table 1. The basic observational properties of the different WFC3/IR survey fields analysed in this work. Column 1 lists the adopted field name, columns 2 
& 3 list the coordinates of the centre of the field and column 4 lists the survey area in arcmin . Columns 5-14 list the global average 5a aperture depths in 
each of the available filters. As discussed in Section 2.3, the chosen aperture diameters are 0.3"for the ACS imaging in the HUDF, HUDF09-1 and HUDF09-2 
fields, and 0.4"for the ACS imaging in all other fields. For the WFC3/IR filters, the chosen aperture diameters are 0.4", 0.4", 0.44", 0.47"and 0.50" for the 
^098, Yl05, Jl25, </l40 & #160 bands respectively. It should be noted that the depths within a given field can vary significantly (e.g. CANDELS GOOD-S 
DEEP) and that there are no coordinates listed for BoRG simply because it consists of a large number of widely-separated pointings. 



Field 


RA(J2000) 


Dec(J2000) 


Area 


#435 


Vfeoe 


«775 


«814 


2850 


^098 


yios 


</l25 


Jl40 


#160 


HUDF 


03:32:38.5 


-27:46:57.0 


4.6 


29.9 


30.4 


30.1 




29.5 




30.1 


29.6 


29.6 


29.6 


HUDF09-1 


03:33:01.4 


-27:41:11.5 


4.4 




29.3 


29.1 




29.0 




28.9 


29.0 




28.6 


HUDF09-2 


03:32:23.4 


-27:42:52.0 


4.5 




29.4 


29.2 


30.0 


29.0 




29.0 


29.1 




28.7 


CANDELS GS-WIDE 


03:32:38.5 


-27:53:36.5 


35.1 


28.2 


28.6 


28.0 




27.8 




27.3 


27.5 




27.0 


ERS 


03:33:05.5 


-27:51:21.6 


38.4 


28.2 


28.6 


28.0 




27.8 


27.6 




27.8 




27.4 


CANDELS GS-DEEP 


03:32:29.8 


-27:47:43.0 


64.6 


28.2 


28.6 


28.0 




27.8 




28.3 


28.1 




27.7 


CANDELS UDS 


02:17:25.7 


-05:12:04.6 


144.5 




27.8 




27.7 








27.2 




27.1 


BoRG 






180.4 




27.7 








28.2 




27.2 




26.2 



both parallel fields (drizzled onto a 30mas pixel scale) and have also 
used our own reduction of the ACS imaging covering HUDF09-2 
originally obtained as part of the HUDF05 campaign (GO- 10632, 
P.I. Stiavelli). For the HUDF and HUDF09-1 we have made use of 
the publicly-available reductions of the original HUDF ACS imag- 
ing (Beckwith et al. 2006) and the HUDF05 ACS imaging respec- 
tively. Finally, we have also made use of the new ultra-deep is 14 
ACS imaging (128 orbits) obtained as parallel observations dur- 
ing the UDF12 campaign which provides a ~ 60% overlap with 
HUDF09-2. 

2.1.3 GOODS-S 

In addition to the HUDF and parallel fields, we have made ex- 
tensive use of the publicly-available WFC3/IR imaging of the 
GOODS-S field provided by the CANDELS survey (Grogin et al. 
2011; Koekemoer et al. 2011). The CANDELS WFC3/IR data in 
GOODS-S covers a total of 24 WFC3/IR pointings which are di- 
vided into DEEP and WIDE sub-regions. The DEEP sub-region 
consists of 15 WFC3/IR pointings each consisting of 3 orbits of 
Y105 imaging and 5 orbits in both J125 and #i6o- The WIDE sub- 
region consists of 9 WFC3/IR pointings each consisting of a single 
orbit of integration in the Y105, J125 and #160 filters. In addition 
to the CANDELS imaging we have also analysed the ERS data in 
GOODS-S which consists of 10 WFC3/IR pointings, each of which 
were observed for 2 orbits in the Y098, J125 and #160 filters (Wind- 
horst etal. 2011). 

The optical data we have employed in GOODS-S is the 
publicly-available v2.0 reduction of the original GOOD-S ACS 
imaging in the B435, Vtso6, *775 & ^850 filters (Giavalisco et al. 
2004). The total area of the overlapping WFC3/IR+ACS cover- 
age we have analysed in GOODS-S is 138 arcmin 2 (excluding the 
HUDF and parallel fields). A summary of the available filters and 
depths is provided in Table 1. 

2.1.4 The UKIDSS Ultra Deep Survey Field 

In order to significantly increase the total areal coverage of our 
dataset, crucial for constraining the bright end of the 2 > 7 galaxy 
luminosity function, we have also analysed the publicly-available 
CANDELS WFC3/IR+ACS imaging in the UKIDSS Ultra Deep 
Survey (UDS; Lawrence et al. 2007). The full CANDELS dataset 
in the UDS consists of 44 WFC3/IR pointings, each featuring 4/3 
of an orbit of i/ieo imaging and 2/3 of an orbit in J125. Along 



with the primary WFC3/IR observations, the CANDELS campaign 
in the UDS also obtained ACS parallels with each pixel typically 
receiving three orbits of integration in isi4 and one orbit in V606- 
However, due to the focal plane separation of the WFC3/IR and 
ACS cameras, only 32/44 of the WFC3/IR pointings are fully cov- 
ered by the parallel ACS imaging. It is from these 32 WFC3/IR 
pointings, covering a total area of ~ 150 arcmin 2 , that we have 
selected our sample of high-redshift candidates in the UDS. 

In addition to the primary HST imaging data, in refining our 
high-redshift sample in the UDS we have also made extensive use 
of a variety of ground-based datasets. Amongst these the three 
most important are new ultra-deep 2' — band imaging obtained with 
Suprime-Cam on Subaru which reaches a depth of z' — 26.5 (5a; 
1.8"diameter apertures), new ultra-deep VLT+HAWK-I Y— band 
observations of the UDS CANDELS region obtained as part of 
the HUGS programme (P.I. A. Fontana) which reach a depth of 
Y — 26.5 (5a; 1.25"diameter apertures) and the latest DR10 re- 
lease of the UDS if— band imaging which reaches a depth of K — 
25.1 (5a; 1.8"diameter apertures). Moreover, in order to clean the 
sample of low-redshift interlopers, we have utilised a stack of the 
BVRi' Subaru imaging of the UDS described in Furusawa et al. 
(2008), which reaches a depth of > 29 (2a; 1.8"diameter aper- 
tures). Finally, in order to further refine our photometric redshift 
solutions we have exploited narrow- band (NB921, \ c = 9210A) 
imaging of the UDS (Sobral et al. 201 1) which reaches a depth of 
2:921 = 26.0 (5a; 1.8"diameter apertures). 

2.1.5 BoRG 

In order to increase our ability to constrain the bright-end of the 
2 = 8 luminosity function we have performed our own reduction 
and analysis (Bowler et al. 2013, in preparation) of the data taken 
by the Brightest of the Reionizing Galaxies survey (BoRG; Trend et 
al. 2011, 2012). BoRG is a HST pure-parallel programme, consist- 
ing of imaging in four filters from WFC3, designed to detect 2 ~ 8 
Lyman-break galaxies as Y 98 drop-outs. Details of the BoRG ob- 
servation strategy can be found in Trenti et. al. (2011), however 
briefly, pure-parallel observations were obtained at multiple sight- 
lines in the Yogs, J125, Him filters from WFC3/IR and one or both 
of Vf306 and Vt>oo with WFC3/UVIS. The exposure times are cho- 
sen to allow the detection of Yfras drop-out galaxies at z ~ 8 based 
on a large Y ms — J125 colour and relatively flat rest-frame spectral 
slope inferred from the J125 — Hieo colour. 

At the time of writing the complete BoRG dataset consists of 
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69 independent fields with a variety of different exposure times. 
However, to homogenize the dataset somewhat, we have restricted 
our analysis to the 41 BoRG fields with 5<r detection limits in the 
range J125 > 27.2 — 27.9 (0.44"diameter aperture), and which 
lie at high galactic latitude. Based on these 41 fields (total area 
180arcmin 2 ) an initial candidate list was constructed using the 
colour-cuts employed by Bradley et al. (2012), but based on our 
own photometry and depth analysis (see below). As with all the 
other survey fields analysed here, these candidates were then anal- 
ysed with our photometric redshift code before being included in 
the luminosity function analysis. To mitigate the increased pho- 
tometric redshift uncertainties introduced by the small number of 
available filters, in our final analysis of the z = 8 luminosity func- 
tion we include only those candidates selected from BoRG brighter 
than M1500 < -20.5. 

2.2 Catalogue production 

Given that the primary focus of this paper is the evolution of the 
galaxy luminosity function at z > 6.5, by which point the Lyman- 
break has been redshifted to an observed wavelength X bs > 
0.9/im, high-redshift galaxy candidate selection was performed ex- 
clusively in the near-IR using the WFC3/IR imaging available in 
each field. 

In order to provide a master catalogue which was as com- 
plete as possible, objects were initially selected using each individ- 
ual near-IR image and from every possible wavelength-contiguous 
stack of near-IR images. For example, based on the UDF12 data, 
four object catalogues were generated with SEXTRACTOR v2.8.6 
(Bertin & Arnouts 1996) using the individual Y105, J125, Jito and 
-H160 images for object detection and a futher six catalogues were 
generated using the following near-IR stacks as the detection im- 
age: Y105 + J125 + Juo + Hieo, Y105 + J125 + Jim, iios + 
J12S + Jim, J125 + J140 + Hieo & ^140 + Hl60- 

From this initial set of 10 object catalogues, a master cata- 
logue was constructed containing every unique object which was 
detected at > 5er significance in any of the detection images. For 
those objects which were present in multiple catalogues, the posi- 
tional information and photometry based on the highest signal-to- 
noise detection was propagated to the master catalogue. Although 
this specific example is within the context of the UDF12 dataset in 
the HUDF, the general selection process was identical in each field, 
notwithstanding differences enforced by the number of available 
filters. 



2.3 Photometry 

When selecting samples of high-redshift galaxies, the choice of 
photometric apertures is inevitably a balance between optimizing 
depth and ensuring that it is possible to derive well-defined and sta- 
ble aperture corrections. The photometry adopted in this study is 
all based on circular apertures, where the choice of aperture in each 
band is tuned to enclose > 70% of the flux of the filter-specific 
point-spread function (PSF). For those fields with data drizzled 
onto a 0.03"/pix grid, the photometry is based on 0.3"diameter 
apertures in the optical ACS bands and 0.40", 0.44", 0.47"and 
0.50"diameter apertures in the Vios, J125, J140 and Hieo bands. 
For those datasets where the data is drizzled onto a 0.06"/pix grid, 
we have adopted 0.40" diameter apertures in the optical ACS and 
iogsAios WFC3/IR bands, together with 0.44" and 0.50" diameter 
apertures in the J125 and f/160 bands. For the purposes of the sub- 
sequent photometric redshift analysis (see below), the photometry 




Figure 1. Examples of our photometric redshift analysis for three objects in 
the HUDF at redshifts z = 6.9, z = 8.1 and z = 8.8 (top to bottom). In 
each plot the blue curve shows the best-fitting high-redshift solution, while 
the orange curve shows the best-fitting alternative low-redshift solution. In 
each case the inset shows the photometric redshift probability density func- 
tion, p(z), which is incorporated into our luminosity-function analysis. The 
upper limits at 3.6/im and 4.5/^m have been derived via our own decon- 
fusion analysis of the ultra-deep IRAC imaging obtained by Labbe et al. 
(2012), using the technique described in McLure et al. (2011). All upper 
limits are lcr. 



in each band was corrected to a constant enclosed flux level (typi- 
cally 82% — 84%) using the curve of growth of the observed PSFs 
in each band/field. 



2.4 Depth analysis 

A crucial element of reliable high-redshift candidate selection is the 
derivation of accurate information regarding the photometric depth 
of each available image. This information is clearly vital for the 
reliable exclusion of low-redshift contaminants, but is also required 
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Figure 2. The photometric redshift distribution of the full sample of high- 
redshift galaxy candidates (robust+insecure) derived from our analysis of 
the eight survey fields listed in Table 1 . The full sample consists of N = 
576 high-redshift candidates selected from a total area of 477 arcmin 2 . 

in order to provide the robust flux measurement errors necessary for 
accurate photometric redshift results. 

For each survey field in turn, accurate maps of which pix- 
els contained significant object flux were constructed by stacking 
SEXTRACTOR segmentation maps for each available filter. Before 
stacking, the individual segmentation maps were dilated in order 
to better capture the extended wings of bright low-redshift objects. 
Based on these image maps, a large number of photometric aper- 
tures (typically 10-200 thousand) were located amongst the frac- 
tion of each image that had been determined to be dominated by 
'sky'. A robust estimator was then used to measure the sigma of 
the distribution of the fluxes contained with these sky apertures and 
thereby determine the empirical depth of the image for a given aper- 
ture diameter. This aperture-to-aperture r.m.s. depth measurement 
is a robust method of determining the actual significance of any 
aperture flux measurement and captures the true noise properties 
of a given image, which are typically underestimated by a simple 
measurement of the pixel r.m.s. Once the true depth of each im- 
age was determined in this fashion, the corresponding weight maps 
were scaled to match the empirical depth measurement. By includ- 
ing the scaled weight maps in the dual-mode catalogue production 
process, it was thereby possible to provide accurate, position de- 
pendent, flux measurement errors for each source. 

2.5 High-redshift galaxy selection 

In order to study the z > 7 luminosity function it is necessary 
to derive reliable catalogues of z > 6.5 galaxy candidates. Con- 
sequently, the master object catalogues for each survey field were 
initially cleaned of low-redshift contaminants by insisting that each 
object remained undetected at the 2a level in each filter short-ward 
of the zs&o waveband. In addition, to exclude the small number of 
objects with ~ 1.5 — 2a detections in multiple blue optical filters, 
each high-redshift candidate was also required to be undetected 



at the 2a level in an inverse-variance weighted stack of all filters 
short- ward of zsso- 



2.5.1 Photometric redshift analysis 

After the initial exclusion of low-redshift contaminants, the mas- 
ter catalogues for each field were processed using our photomet- 
ric redshift code. For a full description of the photometric red- 
shift code the reader is referred to McLure et al. (2011), but we 
provide a brief outline here for completeness. The photometry for 
each galaxy is fitted with a range of galaxy templates, either em- 
pirical spectra or evolving synthetic galaxy-evolution models, with 
the best-fitting galaxy parameters determined via ^-minimization. 
To ensure a proper treatment of the photometric uncertainties, the 
model fitting is performed in flux-wavelength space, rather than 
magnitude-wavelength space, and negative fluxes are included. The 
code fits a wide range of reddening, based on the Calzetti et al. 
(2000) dust attenuation law, and accounts for IGM absorption ac- 
cording to the Madau (1995) prescription. If necessary, the code 
will also fit each high-redshift candidate with the addition of Lya 
emission based on an input prescription for the allowable range of 
rest-frame equivalent widths (EWo). The best-fitting photometric 
redshifts and redshift probability density functions adopted here are 
based on Bruzual & Chariot (2003) stellar population models with 
metallicities of 0.2Zq or Zq, although the exact choice of stel- 
lar population model has little impact on the derived p(z) in most 
cases. Three examples of the results of our photometric-redshift 
analysis are shown in Fig. 1; the three objects illustrated have de- 
rived redshifts of z ~ 6.9, z ~ 8.1 and z — 8.8 (see Table Al). 

Based on the photometric-redshift results, objects were ex- 
cluded if it was impossible to obtain a statistically acceptable so- 
lution at z > 6.5 (typically Xbest — 20), or if the photometric 
redshift probability density function indicated a very low probabil- 
ity that the object is at z > 6.5 (i.e. J^Z^l p(z)dz < 0.05). 

At this point, the catalogues were visually expected in order to 
remove spurious contaminants such as artefacts, diffraction spikes 
and over-deblended low-z objects. Following this final cleaning 
step, and following the procedure of McLure et al. (2011), each 
object was classified as robust or insecure depending on whether 
or not the best-fitting low-redshift solution could be ruled out at 
the 2a confidence level. The objects classified as robust constitute 
our most reliable high-redshift galaxy sample and are presented in 
the Appendix for comparison with the results of other similar stud- 
ies and potential spectrocopic follow-up observations. However, it 
should be noted that all high-redshift candidates which survived the 
full selection process, whether robust or insecure, were included in 
the luminosity function analysis (see below). 



2.5.2 Dwarf star exclusion 

As has been widely discussed in the literature (e.g. Dunlop 2012; 
Bowler et al. 2012), ultra-cool galactic dwarf stars (M, L, & T 
dwarfs) are a potential source of low-redshift contamination of the 
bright-end of the galaxy luminosity function at z > 6. Within the 
context of the current study, the most serious source of concern is 
potential contamination of the z — 7 galaxy luminosity function 
by T-dwarf stars. To deal with this problem the photometry of all 
bright z ~ 7 LBG candidates was analysed using a library of em- 
pirical optical-to-near-IR dwarf star spectra spanning the spectral 
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range MO to T9, taken from the SpeX archiv^j All bright candi- 
dates which returned a high-quality fit with a dwarf star template 
and had a measured WFC3/IR half-light radius which was consis- 
tent with being spatially unresolved were removed from the final 
sample. 



Table 2. The results of our SWML determination of the z = 7 and 2 = 8 
galaxy luminosity functions. Columns 1 and 3 list the luminosity function 
bins adopted at z = 7 and z = 8 respectively (all bins are 0.5 magnitudes 
wide). Columns 2 and 4 list the individual values of tj)/, and their corre- 
sponding uncertainties. 



2.6 Final galaxy sample 

The final high-redshift galaxy sample utilised in this analysis com- 
prises a total of N = 576 galaxy candidates (robust+insecure) se- 
lected from a total survey area of 477arcmin 2 . The photometric 
redshift distribution of this sample is shown in Fig. 2. The indi- 
vidual redshift probability distributions, p(z), of all these galaxies 
were used in the LF determination. In Table Al in the Appendix 
we provide a catalogue of the 100 most robust 2 ~ 6.5 — 12.0 
galaxies uncovered within the HUDF itself (i.e. from the UDF12 
data), with positions, photometric redshifts (including uncertain- 
ties), total absolute magnitudes Afisoo, apparent Hxeo magnitudes 
(0.5"aperture), and cross-referencing to previous studies as appro- 
priate. The ACS+WFC3/IR photometry for these same objects is 
provided in Table A2. Tables for the robust objects in the other 7 
survey fields analysed here will be provided in the online version 
of the published paper. 

We note that the advance at faint magnitudes offered by the 
UDF12 dataset is clear, because 50 of the 100 robust 2 > 6.5 
HUDF galaxies tabulated in Table Al have not been reported in 
any previous study. Of these additional 50, colour-colour selection 
(as described by Schenker et al. 2013) finds 23 sources. The re- 
maining 27 galaxy candidates are only revealed by our photometric 
redshift analysis which exploits all of the 4-band WFC3/IR imag- 
ing. Note that, for simplicity, we have decided to report only -H160 
magnitudes in Table Al, and so some robust objects may appear to 
be surprisingly faint in i^ieo because they are better detected in the 
shorter-wavelength WFC3/IR filters. 

2.7 Simulations 

The final result of the object selection process is a catalogue of can- 
didate high-redshift galaxies, each with an associated photometric 
redshift probability density function. However, in order to accu- 
rately derived the galaxy luminosity function it is then necessary to 
employ detailed simulations to map between the derived and intrin- 
sic properties of each candidate. 

Following the methodology of McLure et al. (2009), we adopt 
a parametric model of the evolving high-redshift galaxy luminos- 
ity function in order to generate a realistic synthetic population 
of high-redshift galaxies. The model adopted here is a Schechter 
function with the three parameters ((j>* , M^y, a) evolving linearly 
with redshift, changing from (9.8 x 10" 4 Mpc" 3 , -20.7, -1.65) 
at z = 4.5 to (2.3 x lO^Mpc" 3 , -20.1, -2.1) at 2 = 9.0. Al- 
though simple, this parameterization successfully reproduces the 
observed evolution of the UV-selected galaxy luminosity function 
within the redshift interval 5 < 2 < 8 (e.g. McLure et al. 2009, 
2010; Bouwens et al. 2007, 201 1; Bradley et al. 2012). 

For each survey field, the evolving luminosity function model 
is used to populate an input apparent magnitude - redshift grid 
(wii50u — z, where misoo is the apparent mag nitude at 1500A), 
which is divided into cells of width Sz — 0.05 and 5mi5oo = 0.1. 
For each simulated object, synthetic photometry is generated using 

B http://pono.ucsd.edu/ adam/browndwarfs/spexprism/ 



z = 7 z=8 





</> fe /mag 1 Mpc 3 


M1500 


</> fe /mag J Mpc 3 


-21.00 
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-18.75 
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-18.00 


0.0042±0.0009 


-18.25 


0.0018±0.0006 


-17.50 


0.0079±0.0019 


-17.75 


0.0028±0.0008 


-17.00 


0.011±0.0025 


-17.25 
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Bruzual & Chariot (2003) galaxy templates with a range of redden- 
ing designed to produce a population of objects whose distribution 
of UV slopes is centred on /3 = —2 with a small dispersion (c.f. 
Dunlop et al. 2012, 2013; Rogers, McLure & Dunlop 2012). The 
synthetic galaxies were then injected into the real data, using the 
empirical measurement of the WFC3/IR and ACS PSF appropri- 
ate for each individual filter. It should be noted that simulations 
were also conducted in which the synthetic galaxies were modelled 
as spatially resolved, with half-light radii drawn from the distri- 
bution measured for 2 = 7 — 8 LBGs by Oesch et al. (2010b). 
However, these simulations were not adopted because they were 
found to provide results virtually identical to the injection of PSFs 
and in the regime where the simulation results are most crucial for 
the luminosity function determination (i.e. Afisoo > —19) several 
studies indicate that the UV-selected galaxy population is virtually 
unresolved at the resolution of WFC3/IR (i.e. Oesch et al. 2010b; 
Grazian et al. 2012; Ono et al. 2013). 

The simulated galaxy population was then selected and pro- 
cessed through our photometric redshift analysis in an identical 
fashion to the real sample of high-redshift galaxy candidates. The 
resulting p(z) distribution for each synthetic galaxy was then used 
to populate an output misoo — z grid. Under the assumption 
that the simulation provides a reasonable description of the actual 
high-redshift galaxy population, the ratio of the output and input 
W.1500 — z grids provides a mapping between the observed and in- 
trinsic properties of the high-redshift galaxy population, which au- 
tomatically accounts for selection efficiency, photometric redshift 
errors and flux boosting. 



3 LUMINOSITY FUNCTION ESTIMATION 

For each survey field, the p(z) distributions for each high-redshift 
galaxy candidate were used to populate an observed mi 5 oo — z 
plane, which was then corrected using the results of the correspond- 
ing simulation to provide our best estimate of the distribution of the 
observed high-redshift population on the intrinsic misoo — 2 plane. 
It is this information that is then used to estimate the galaxy lumi- 
nosity function using two different techniques. 

The primary method is an implementation of the non- 
parametric step- wise maximum-likelihood (SWML) method of Ef- 
stathiou, Ellis & Peterson (1988). It is the results of this method 
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Figure 3. The left-hand panel shows our new determination of the UV-selected galaxy luminosity function at redshift z = 7. The data-points have been derived 
using a combination of our photometric redshift analysis and an implementation of the step-wise maximum likelihood (SWML) method of Efstathiou, Ellis & 
Peterson (1988); see text for details. The thick solid red line shows the best-fitting Schechter function derived using the parameteric STY maximum-likelihood 
technique applied simultaneously to the first seven survey fields listed in Table 1. For comparison, we also show a straight-forward x 2— fit to the binned 
SWML data (thick black line). The right-hand plot shows the same information for our new determination of the z = 8 luminosity function. To derive the 
2 = 8 luminosity function we have also incorporated the information derived from our own reduction and analysis of the BoRG dataset. 



which provide our basic determination of 4>(Mi500, z), without re- 
lying on the assumption that the luminosity function obeys a par- 
ticular functional form. However, in order to compare with previ- 
ous results and to study the redshift evolution of the galaxy lumi- 
nosity function, it is also desirable to derive parametric fits to the 
galaxy luminosity function. Therefore, in order to derive Schechter- 
function fits to the galaxy luminosity function we have also imple- 
mented a version of the parametric maximum-likelihood technique 
of Sandage, Tammann & Yahil (1979); hereafter STY. 

When defining the likelihood, both techniques rely on the as- 
sumption that, at a given redshift, the probability of observing a 
galaxy of a given luminosity can be defined as follows: 



Pi oc 



<t>{U) 



£«" 4>(L)dL 



(1) 



where Ln m is the limiting luminosity of the survey. Ideally this 
should be implemented in the situation where each galaxy has a 
unique spectroscopic redshift and luminosity. In the absence of 
this information, our implementation does the next best thing and 
adopts the normalized probability density function for each high- 
redshift candidate, with the absolute UV magnitude (Misqo), cal- 
culated using a top-hat filter at 1500A in the rest-frame of the best- 
fitting SED template, re-calculated at each step within the p(z). The 
overall best-fit is determined by maximizing the following likeli- 
hood: 



(2) 



where the outer product symbol indicates that the maximum likeli- 
hood is calculated over j separate survey fields, each with its own 
limiting luminosity. 

One of the key strengths of both techniques is that they take 
no account of the absolute number density of objects and should 
therefore be relatively insensitive to the effects of cosmic struc- 
ture. However, as a result, it is necessary to determine the overall 



normalization of the resulting galaxy luminosity-function estimates 
independently. In each case, we have derived the luminosity func- 
tion normalization by requiring that the maximum-likelihood lumi- 
nosity function estimates reproduce the cumulative number counts 
of observed galaxies within the appropriate redshift and absolute- 
magnitude intervals. 



4 THE GALAXY LUMINOSITY FUNCTION RESULTS 

In Fig. 3 we show our new determinations of the UV-selected 
galaxy luminosity functions at z — 7 and z = 8. In each panel 
the data-points have been derived using the SWML technique, and 
include error estimates which have been derived via boot-strap re- 
sampling of the underlying high-redshift galaxy sample. In Table 2 
we provide the individual, step-by-step, SWML determinations of 
the z — 7 and z — 8 luminosity functions and their corresponding 
uncertainties. 

In each panel of Fig. 3 the thick red line shows our best-fitting 
Schechter function as derived via our implementation of the STY 
maximum-likelihood technique, which is constrained via a simul- 
taneous fit to all relevant survey fields (as listed in Table 1). The 
confidence intervals on the derived faint-end slope and characteris- 
tic magnitude from our STY fits at z = 7 and z — 8 are shown in 
Fig. 4 and the best-fitting Schechter-function parameters are listed 
in Table 3. 

As can be seen from the results presented in Fig. 3, the ad- 
ditional dynamic range in luminosity provided by combining the 
new UDF12 dataset with the wider area GOODS-S, CANDELS- 
UDS and BoRG datasets has allowed an accurate determination 
of the 2 = 7 and 2 = 8 luminosity functions spanning a factor 
> 50 in UV luminosity. In particular, the new UDF12 dataset has 
allowed us to constrain the form of the luminosity function as faint 
as L < 0.1L* for the first time. As can be seen from the results 
presented in Table 3, our new analysis has confirmed that the lumi- 
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Figure 4. The confidence intervals for the faint-end slope (a) and char- 
acteristic absolute magnitude (M* B00 ) derived from our STY maximum- 
likelihood fits to the galaxy luminosity function at z = 7 (top) and z = 8 
(bottom). In each plot the lcr and 2cr confidence intervals are illustrated by 
the dark and light-blue shaded areas respectively. It should be noted that the 
confidence intervals plotted are based on likelihood ratios and correspond 
to Ax 2 = 1 & 4 from the best overall fit (shown by the filled circle in each 
case). These specific confidence intervals have been chosen in order that the 
corresponding one parameter uncertainties can be calculated by projecting 
the contours onto the relevant axis. 

nosity function remains extremely steep at Afisoo > —18 at both 
z — 7 and z = 8, being consistent with a = —2.0 in both cases. 
Moreover, our new analysis confirms and strengthens previous re- 
sults which suggested that there is little evolution in M* 500 between 
z = 7 and z — 8. In the next section we briefly compare our re- 
sults with other relevant results in the literature, before proceeding 
to consider the evolution of the luminosity function from z ~ 6 to 
z ~ 10 and its implications for cosmic reionization. 



4.1 Comparison with previous results 

Although an exhaustive comparison with previous high-redshift lu- 
minosity function work in the literature is beyond the scope of this 
paper, it is instructive to compare the results derived here with 
those of other recent studies. Consequently, we will concentrate 
on a comparison between our new results and those of other HST 
studies which have attempted to fit Schechter-f unction parameters 
and their corresponding uncertainties. The details of the Schechter- 
function fits derived by the various different studies are provided in 
Table 3. 

Given that it also includes the new UDF12 data, it is obvi- 
ously of interest to compare our results to those of the companion 
drop-out analysis of the z — 7 and z — 8 luminosity functions 
performed by our team (Schenker et al. 2013). It can be seen from 
Table 3 that at z = 7 our results are in good agreement, with the 
extra luminosity leverage provided by UDF12 leading both stud- 
ies to conclude that the faint-end slope is a = —1.9, with a small 
uncertainty. The overall uncertainties on the Schechter-function pa- 
rameters derived in this work are slightly tighter than in Schenker 
et al. (2013), which is expected given that we have analysed a larger 
survey area. At z = 8 our results are again consistent with those 
published in Schenker et al. (2013), particularly in term of the faint- 
end slope, although the value of M*soo derived by Schenker et al. 
is somewhat brighter than found here. 

Before the start of the UDF12 imaging campaign, the most 
comprehensive study of the z — 7 and z — 8 luminosity func- 
tions was performed by Bouwens et al. (2011a), who combined a 
dropout analysis of the WFC3/IR imaging in the HUDF09 and ERS 
(total area 53 arcmin 2 ) with constraints provided by various wider- 
area datasets. It can be seen from Table 3 that, in terms of derived 
luminosity function parameters, there is actually very good agree- 
ment between the new results derived here and those of Bouwens 
et al. (2011a). At some level this may be slightly fortuitous, given 
that the agreement between the Schechter function parameters is 
better than that between our respective binned SWML results and, 
as noted by Bouwens et al. (2011a), the small area of their study 
meant that their z — 8 faint-end slope could have been as shallow 
as a — —1.67 ± 0.40 depending on the inclusion/exclusion of two 
bright candidates in HUDF09-2. However, irrespective of this, it 
is clear that the fundamental advantage of our new analysis is the 
ability to better constrain the faint-end slope at both redshifts, using 
the deeper UDF12 imaging. 

Much of the recent work in the literature has been focused on 
trying to improve our knowledge of the 2 = 8 luminosity function. 
Within this context, two recent studies by Oesch et al. (2012) and 
Bradley et al. (2012) have investigated the form of the z = 8 lu- 
minosity function by combining the faint-end results of Bouwens 
et al. (2011a), with improved constraints at the bright end. In the 
case of Bradley et al. (2012) the bright-end information is provided 
by their drop-out analysis of the BoRG WFC3/IR parallel observa- 
tions, whereas in Oesch et al. (2012) the bright-end constraints are 
provided by a drop-out analysis of the CANDELS DEEP+WIDE 
imaging in GOODS-S (total area 95 arcmin 2 ). It can be seen from 
the results presented in Table 3 that, in terms of derived luminosity- 
function parameters, there is very good agreement between the new 
results derived here and those of Oesch et al. (2012) and Bradley et 
al. (2012), with all studies seemingly converging on a steep faint- 
end slope of a — —2.0 and M*soo — —20.1. 

Overall, the comparison between derived luminosity function 
parameters shown in Table 3 is therefore highly encouraging, es- 
pecially given the different datasets, reductions and analysis tech- 
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Table 3. A comparison between the Schechter-function parameters for the 
z = 7 and 2 = 8 galaxy luminosity function derived by various recent HST 
studies. The first column lists the name of the study and columns two to 
four list the Schechter-function parameters and their quoted uncertainties. 
The units of <f>* are Mpc~ li . 



C\2 
I 



Study 


-^1500 


log(0*) 


a 


z = 7 








This work 


-i9.ootg:l2 


-2 qfi+ - 18 
z.ao_ 23 




Schenker et al. (2013) 


-20.14±g;»S 


°' ly -0.24 


-i.»c_ 17 


Bouwens et al. (201 la) 


-20.14+°;" 


, n7 +0.26 
J - u '-0.26 


„ Ql+0.21 
z - ul -0.21 



This work 


—20 1 2+ 0,37 
zu.iz_ Q 4g 


„ oc+0.28 
o - OJ -0.47 


-2.02 


Schenker et al. (2013) 


-20 44+ 47 


q 5n +0.35 


-1.94 


Bouwens et al. (2011a) 


-20.iol°;g 


q r.q+0.43 

°- zo -0.43 


-1.91 


Oesch et al. (2012) 


-M.M+g;" 


O.OU_q 4 g 


-2.06 


Bradley et al. (2012) 


-20.26lg;» 


q q 7 +0.26 

• 3 -°'-0.29 


-1.98 



+0.22 
-0.23 



+0.21 

-0.24 



+0.32 
-0.32 



+0.35 



+0.23 
-0.22 



niques adopted by the various different studies. However, due to the 
fact that the current study (together with Schenker et al. 2013) ex- 
ploits the deeper imaging provided by UDF12 and, uniquely, incor- 
porates WFC3/IR imaging covering a wider area than all previous 
studies (including GOODS-S, CANDELS-UDS and BoRG) we are 
confident that the luminosity-function determination provided here 
is the most accurate currently available at these redshifts. 

4.2 The evolution of the luminosity function 

Several previous studies have concluded that the evolution of the 
galaxy luminosity function over the redshift range 5 < z < 7 can 
be well described as pure luminosity evolution (e.g. Bouwens et 
al. 2007; McLure et al. 2009; Bouwens et al. 2011a). It is clearly 
of some interest to investigate whether or not the evolution of the 
luminosity function from 2 = 7 — 8 remains consistent with this 
apparently simple picture. 

Some insight into this question can be gained by examining 
the confidence intervals on the faint-end slope and characteristic 
magnitude at z = 7 and 2 = 8 shown in Fig. 4. It can imme- 
diately be seen from Fig. 4 that our new analysis provides little 
evidence for a significant change in Mi 50 o or a over the redshift 
inteval 2 = 7 — 8 and in fact, the best-fitting Schechter-function pa- 
rameters (see Table 3) suggest that the dominate change is a factor 
of ~ 2.5 drop in <j>* between 2 = 7 and 2 = 8. We note here that 
Bouwens et al. (201 la) also commented that some of the 2 = 7—8 
evolution may be explained by a change in <f>*, but concluded that 
the uncertainties were too large to be confident. Although our im- 
proved determinations of the 2 = 7 and 2 = 8 luminosity functions 
strengthen the suggestion that 0* is changing within the redshift 
range 7 < 2 < 8, the results presented in Fig. 5 indicate that the 
available data are still insufficient to rule out pure luminosity evo- 
lution. 

Given that pure luminosity evolution provides such a good de- 
scription of the luminosity function evolution over the redshift in- 
terval 4 < z < 7 (Bouwens et al 2011a), in Fig. 5 we explore 
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Figure 5. An illustration that pure luminosity evolution can provide an ac- 
ceptable fit to the observed evolution of the galaxy luminosity function over 
the redshift interval 6 < z < 8. The blue and red data-points show the 
S WML determination of the galaxy luminosity function at z = 7 and 2 = 8 
from this work, while the green data-points show the determination of the 
2 = 6 luminosity function from McLure et al. (2009). The corresponding 
curves show the results of fitting the binned LF data with Schechter func- 
tions where the faint-end slope and overall normalization have been held 
constant at representative values (a = —1.9 and (/>* = 0.00085 Mpc~ 3 
respectively) but the characteristic magnitude (.M* SQ0 ) has been allowed to 



float. In this scenario Mf_ 500 evolves by ~ 0.3 magnitudes per Az 



interval, changing from Mf 



-20.3 at 2 = 6 to Mf 



1 

-19.7 



at 2 = 8. It can be seen that this simple parameterization is capable of 
satisfactorily reproducing the observed data. 

whether a simple luminosity evolution parameterization can con- 
tinue to provide an adequate description of the observed evolution 
at 2 > 6. To investigate this issue we simply fit Schechter functions 
to the binned SWML luminosity-function data at 2 = 6,7&8, 
allowing Mxsoo t0 A° a t as a f ree parameter, but keeping cf>* and 
a fixed at representative values (a = —1.9 and (j>* = 0.00085 
Mpc~ 3 respectively). In this simplied scenario we find that Mf 500 
evolves by ~ 0.3 magnitudes per Az = 1 interval, changing from 
M* 500 ~ -20.3 at 2 = 6 to M* 500 ~ -19.7 at 2 = 8. It is imme- 
diately clear from Fig. 5 that, within the constraints of the current 
data, it is still perfectly possible to reproduce the observed lumi- 
nosity function data in the redshift interval 6 < z < 8 with pure 
luminosity evolution alone. 



4.3 The galaxy luminosity function at z = 9 

In addition to the quadrupling of the available Yios imaging in the 
HUDF, the key advantage provided by the new UDF12 dataset is 
the addition of ultra-deep imaging in the previously unexploited 
J140 filter. The availability of the new J140 imaging provides the 
first real opportunity to constrain the faint end of the 2 ~ 9 lumi- 
nosity function, simply because, in the redshift interval 8.5 < 2 < 
9.5, the J140 and iifieo imaging still provide two filters long-ward 
of the redshifted Lyman break. 

In Fig. 6 we show our SWML determination of the 2 = 9 
galaxy luminosity function, which is derived entirely from the data 
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Table 4. The results of our SWML determination of two luminosity bins 
on the z = 9 galaxy luminosity function. Column 1 lists the adopted lumi- 
nosity function bins (both are 0.5 magnitudes wide) and column 2 lists the 
individual values of <j>^ and their corresponding uncertainties. 



Afl500 


fc /mag iJVTpc 3 


-18.00 


0.0016±0.0007 


-17.50 


0.0021±0.0009 



available in the HUDF itself. For comparison, we also show in Fig. 
6 our SWML determination of the z = 8 luminosity function, 
which is identical that shown in Fig. 3. Although it is clearly not 
sensible to draw strong conclusions from two low signal-to-noise 
luminosity-function bins, the z — 9 data-points shown in Fig. 6 
(and listed in Table 4) immediately suggest that there is no dra- 
matic fall in the volume density of M1500 — —18 galaxies between 
z — 8 and 2 = 9. The solid and dashed lines in Fig. 6 show the 
result of Schechter function fits to the two z = 9 data-points under 
the assumption that the evolution from z — 8 to z = 9 is either 
purely luminosity evolution (solid line) or purely density evolution 
(dashed line). In both cases it is assumed that the luminosity func- 
tion faint-end slope remains unchanged at a = —2.02. 

It can clearly be seen from Fig. 6 that it is impossible to say 
anything meaningful about the form of any evolution of the galaxy 
luminosity function between z — 8 and z = 9. However, the two 
alternative Schechter-function fits do demonstrate that, provided 
the faint-end slope does not change substantially between z — 8 
and z = 9, irrespective of the form of the evolution the resulting 
integrated UV luminosity density is likely to be very similar. 



5 EVOLVING LUMINOSITY DENSITY AND COSMIC 
STAR-FORMATION HISTORY 

Given our improved knowledge of the evolving galaxy luminosity 
function, it is clearly of interest to briefly explore the implications 
for the evolution of the observed UV luminosity density, which has 
a direct impact on determining whether the observed high-redshift 
galaxy population can reionize the Universe. For a thorough re- 
view of the constraints which could be placed on reionization pre- 
UDF12, the reader is referred to Robertson et al. (2010) and Finkel- 
stein et al. (2012). 

In Fig. 7 we provide a new calculation of the evolution of the 
observed UV luminosity density based on the McLure et al. (2009) 
determination of the z = 6 luminosity function and the new lumi- 
nosity function determinations derived in this work at z = 7 — 9. 
In all cases the datapoints show the results of integrating the ap- 
propriate luminosity function down to an absolute magnitude of 
M1500 = —17.7, in order to allow straightforward comparison 
with previous work. In both panels the right-hand axis shows how 
the UV luminosity density converts into a star-formation rate den- 
sity under the assumption of a Salpeter IMF and the conversion be- 
tween UV luminosity and star-formation rate prescribed by Madau, 
Pozzetti & Dickinson (1998). In the top panel we show the evolu- 
tion in UV luminosity density as a function of redshift, whereas in 
the bottom panel we show the evolution as a function of look-back 
time. The upper and lower datapoints at z = 10 show how the 
UV luminosity density changes if the luminosity function evolu- 
tion from z = 8 — 9 continues to z = 10 as either pure luminosity 




-22 -20 -18 

M 1500 

Figure 6. The large purple data-points show our SWML estimate of the 
z ~ 9 luminosity function derived from the UDF12 dataset. For context, 
in grey we also show our estimate of the z = 8 luminosity function. The 
solid and dashed lines show fits to the z = 9 data-points under the assump- 
tion that the luminosity function evolves from z = 8 to z = 9 via pure 
luminosity (solid) or pure density (dashed) evolution respectively (both as- 
sume the faint-end slope remains fixed at the 2 = 8 value). Although it is 
currently impossible to differentiate between them, it is clear that both sce- 
narios would lead to very similar integrated UV luminosity densities (see 
text for a discussion). 

(lower) or pure density evolution (upper) respectively. No dust cor- 
rections have been applied to any of the datapoints. 

It can be seen from Fig. 7 that our data indicate that there is an 
order of magnitude increase in the observed UV luminosity density 
over the 450 Myr period between z = 10 and z = 6. Moreover, 
from the bottom panel of Fig. 7 it can be seen that the increase 
in UV density is very close to linear with cosmic time between 
z = 8 and z = 6. However, our new information at z ~ 9 from 
the UDF12 dataset provides some evidence that the fall-off in UV 
density at z > 8 is steeper than a linear trend with cosmic time, 
particularly if the galaxy luminosity function continues to evolve 
primary via luminosity evolution. We note that in this regard our 
results are consistent with the recent results based on the CLASH 
clusters campaign published by Coe et al. (2013) and Bouwens et 
al. (2013), which are plotted in the bottom panel of Fig. 7. A full 
discussion of the implications of these results within the context of 
cosmic reionization will appear in Robertson et al. (2013, in prepa- 
ration). 



6 CONCLUSIONS 

By combining the extreme near-IR depth provided by the UDF12 
campaign with extensive wider area WFC3/IR imaging data, it has 
been possible to study the high-redshift galaxy luminosity function 
using galaxy samples selected from a total area of ~ 480 arcmin 2 , 
which span a factor of > 50 in luminosity and a factor of ~ 1000 in 
number density. Based on this unique dataset it has been possible to 
determine the most accurate measurement to date of the jz = 7 and 
z — 8 galaxy luminosity functions and the first meaningful con- 
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fitting Schechter functions, our formal constraints on the faint-end 
slope are a = — 1.90jlo'.i5 at 2 = 7 and a - 



z - uz -0.23 



at z 




12.6 12.8 
look— back time/Gyr 



Figure 7. The evolution of the observed UV luminosity density as a func- 
tion of redshift (top) and look-back time (bottom). The data-points at 
6 < z < 9 show the results of integrating the galaxy luminosity function 
to an absolute magnitude limit of Afisoo = —17.7, where the luminos- 
ity function parameters are taken from McLure et al. (2009) at z = 6 and 
from this work for 2 = 7 — 9. The upper and lower data-points at z = 10 
are derived by assuming that the evolution of the luminosity function from 
2 = 8 — 9 continues as either pure luminosity (lower) or pure density (up- 
per) evolution respectively. The dashed line shown in the bottom panel is 
a linear fit to the evolution of the UV luminosity density between z = 6 
and z = 8. In the bottom panel the two pink data-points are the recent 
UV luminosity estimates from Coe et al. (2013), while the red data-point 
is the estimate from Bouwens et al. (2013). In each panel the right-hand 
axis shows how UV luminosity density converts to star-formation density 
assuming a Salpeter IMF and the UV-to-SFR conversion of Madua, Pozzetti 
& Dickinson (1998). 



(ii) The results of our Schechter function fits strengthen previ- 
ous suggestions that the form of the evolution of the luminosity 
function between 2 = 7 and z = 8 is more akin to density evo- 
lution, rather than the apparent luminosity evolution observed at 
redshifts z = 5 — 7. 

(iii) However, even with the extra leverage provided by the 
UDF12 dataset, we conclude that it is not possible to differenti- 
ate between luminosity and density evolution between z = 7 and 
z — 8. In fact, we demonstrate that it is perfectly possible to pro- 
vide an adequate description of the observed luminosity function 
data between z = 6 and z = 8 under the assumption of pure lumi- 
nosity evolution alone. 

(iv) The unique nature of the UDF12 dataset has allowed us to 
place the first meaningful constraints on the faint end of the galaxy 
luminosity function at z = 9. Taken at face value, these initial re- 
sults suggest that, at least at Mi 500 — —18, there is not a dramatic 
fall-off in the volume density of faint galaxies between z = 8 and 
z = 9. 

(v) Based on our determinations of the galaxy luminosity func- 
tion within the redshift interval 6 < z < 9, we briefly explore the 
evolution of the observed UV luminosity density. Our results indi- 
cate that there is an order of magnitude increase in the UV lumi- 
nosity density over the redshift range 6 < 2 < 10 and that between 
z — 8 and z = 6 the UV luminosity density increases linearly 
with cosmic time. However, our new results at z ~ 9, together with 
recent results from the literature, suggest that the fall-off in UV lu- 
minosity density at 2 > 8 is steeper than would be expected for a 
linear trend with cosmic time. 
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straints on the galaxy luminosity function at 2 ~ 9. The principal 
results and conclusions of our study can be summarised as follows: 

(i) The extra depth provided by the UDF12 dataset has allowed 
us to demonstrate that the faint-end slope of the galaxy luminosity 
function at 2 = 7 and 2 = 8 remains extremely steep, down to 
M1500 = —16.75 and M1500 = —17.00 respectively. Based on 
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APPENDIX A: ROBUST HIGH-REDSHIFT CANDIDATES 
IN THE HUDF 

In Table Al we provide the coordinates, photometric redshifts and 
absolute UV magnitudes (A/1500) for our N=100 robust z > 6.5 
galaxy candidates in the HUDF field. In Table A2 we provide the 
z 850, ^105, ^125, J140 an d -H160 photometry for the same sample. 
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Table Al. Candidate z > 6.5 galaxies in the HUDF. Column one lists the candidate names and columns two and three list the coordinates. Columns four and 
five list the best-fitting photometric redshift and the corresponding lcr uncertainty. Column six lists the total absolute UV magnitude, measured using a top-hat 
filter at 150()A in the rest-frame of the best-fitting galaxy SED template. Column seven lists the -f/ieo aperture magnitude (0.5"diameter = 70% enclosed flux 
for the PSF). Column eight gives references to previous discoveries of objects: (1) McLure et al. (2010), (2) Oesch et al. (2010a), (3) Finkelstein et al. (2010), 
(4) Bunker et al. (2010), (5) Yan et al. (2010), (6) Bouwens et al. (2010), (7) Wilkins et al. (2011) (8) Lorenzoni et al. (2011), (9) Bouwens et al. (2011a), 
(10) Bouwens et al. (2011a) potential, (11) McLure et al. (201 1), (12) Bouwens et al. (201 lb) (13) Schenker et al. (2013). 
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-18 


.5 


28.8+J, I 


9,11,13 


UDF12-4122-7232 


03 


:32 


41 


.22 


-27 


:47: 


:23.2 


6.8 


6.5- 


-7.0 


-17 


.7 


29.7±°;| 




UDF12-3894-7456 


03 


:32 


:38 


.94 


-27 


:47: 


45.6 


6.8 


6.4- 


-7.1 


-17 


.3 


29.9l° ; a 3 




UDF12-4290-7174 


03 


:32: 


:42 


.90 


-27 


:47: 


:17.4 


6.8 


6.4- 


-7.2 


-17 


.3 


30.6l°;I 




UDF12-4056-6436 


03 


:32: 


:40 


.56 


-27 


:46: 


:43.6 


6.8 


6.7- 


-7.0 


-18 


.7 


28.71°;! 


1,2,3,4,5,7,9,11,13 


UDF12-4431-6452 


03 


:32: 


:44. 


.31 


-27 


:46: 


:45.2 


6.8 


6.6- 


-7.0 


-18 


.7 


28.81°;! 


1,9,11,13 


UDF12-3958-6565 


03 


:32 


:39 


.58 


-27 


:46: 


:56.5 


6.8 


6.6- 


-7.0 


-18 


.9 


28.41°;! 


1,2,3,5,9,11,13 


UDF12-4037-6560 


03 


:32 


:40 


.37 


-27 


:46: 


56.0 


6.8 


6.5- 


-7.1 


-17 


.5 


30 l +0 - 4 
JU -0.3 


9,13 


UDF12-4019-6190 


03 


:32 


:40 


.19 


-27 


:46: 


:19.0 


6.9 


6.5- 


-7.2 


-17 


.5 


29.8±°;| 


13 


UDF12-4422-6337 


03 


:32 


:44 


.22 


-27 


:46 


:33.7 


6.9 


6.4- 


-7.3 


-17 


.1 


30.6l ;| 




UDF12-4472-6362 


03 


:32: 


:44 


.72 


-27 


:46: 


:36.2 


6.9 


6.5- 


-7.1 


-18 


.3 


28.71°,! 


13 


UDF12-4263-6416 


03 


:32: 


:42 


63 


-27 


:46: 


:41.6 


6.9 


6.5- 


-7.2 


-17 


.3 


30.2l°; 4 3 
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The z = 7 — 9 galaxy luminosity function 



Table Al. Candidate z > 6.5 galaxies in the HUDF continued.. 



Name 


RA(J2000) 


Dec(J2000) 


^phot 




Misoo 




References 


TTI~\C10 A AQ A 

U Dr 1 2-4484-6568 


ftl.OO.yl /I QA 

03.32.44.84 


— 2/. 46. 56. 8 


£. n 
6.9 


6.5- 


7 1 
- /.I 


— 17.9 


oo ,:+0.4 
29 - 6 -0.3 




UUF12-3975-7451 


AO O r\ OA "7C 

03:32:39.75 


—27:47:45.1 


6.9 


A /I 
6.4- 


7 1 

- /.I 


— 18.1 


7A A +0.2 

29 - 4 -0.2 


9,13 


nr\ci o Qnon <ci on 

uur 12-3989-6189 


03:32:39.89 


T7 . /I ^ . 1 O A 

— 27:46:18.9 


6.9 


6.5- 


7 1 

- /.I 


— 18.0 


7Q c+0.2 

29.5_ 02 


1,13 


UUr 12-3729-61 75 


03:32:37.29 


— 27:46:17.5 


6.9 


6.5- 


7 1 

- /.I 


— 18.2 


7Q n+0.1 

28 - 9 -o.i 


10 


UUF 12-3456-6494 


03:32:34.56 


—27:46:49.4 


7.0 


6.6- 


7 Q 


— 17.9 


on 7+0.1 
29 - 2 -0.1 


9,13 


UDr 12-4068-6498 


03:32:40.68 


— 27:46:49.8 


7.0 


O. J- 


7 A 
- / .4 


— 17.9 


nn ft +0.2 

29.6_ 02 


13 


UUr 12-3692-65 16 


03.32.36.92 


— 2/. 46. 51. 6 


7 A 
/.0 


o.o- 


7 A 


1 7 < 

— 1 /.5 


QA Q-f-0.5 

30.3I 3 


1 n 1 q 
10,1 J 


UDr 1 2-407 1-7 347 


03:32:40.71 


n . -i -7.o /i "7 

— 27:47:34./ 


7.0 


6.5- 


7 Q. 
- / .J 


— 17.8 


QA 1+0-4 


9,13 


TTT~\T7MO AM£. OAOO 

UDr 12-4036-8022 


AQ .TT./lfl q /r 

03:32:40.36 


— 27:48:02.2 


7.0 


6.6- 


7 /I 
— / .4 


— 17.3 


30.6_ 5 


3,9,13 


UDM2-3/55-6019 


AQ .QO.Q"7 CC 

03:32:3 /.55 


— 27:46:01.9 


7.1 


o. /- 


7 /I 
- / .4 


— 17.6 


7Q ft+0-2 

29.6I 02 


9,13 


UUr 12-4256-6566 


03:32:42.56 


—27:46:56.6 


7.1 


"7 ft 


7 7 
- 1 .1 


—20.3 


1£ Q + 0.0 

26 - 9 -o.o 


1,2,3,4,5,7,9,11,13 


T ir\Dl O /I 1 AC "7 1 C/. 

U Dr 1 2-4 1 05- / 1 5 6 


AQ .TT./I 1 AC 

03:32:41.05 


— 27:47:15.6 


7.1 


6.5- 


7 Q 


— 19.0 


7Q /I+0.1 

28 - 4 -o.i 


1,2,3,5,13 


T TI~\C1 1 Q O CQ TCI n 

UUr 12-3853-/519 


AQ .n.TO CQ 

03:32:38.53 


T7 . A -7, c 1 A 

—27:47:51.9 


7.1 


o.y- 


7 3 


— 18.0 


29 ' 7 -0.2 


2,4,13 


T TTVD1 O TOiC /CC^/C 

UUr 12-3825-6566 


AQ .QO.QO TC 

03:32:38.25 


—27:46:56.6 


7.1 


o.o- 


7 A 
— / .0 


— 17.5 


o_a o-|-0 . 3 

29.8_ 02 




UDr 12-3836-61 19 


AQ.QO-QO QA 

03.32.38.36 


— 2/. 46.1 1.9 


7.1 


0.5- 


7 "3 
— /.J 


— 18.8 


28.6_ A 


i a c n 1 1 
3,4,5,9,1 1 


UDr 12-3709-6441 


AQ . QO.Q"7 A A 

03:32:37.09 


n, A£. A A 1 

— 27:46:44.1 


7.2 


6.6- 


7 A 

- / .0 


— 17.0 


QA 7 + 0-4 

30 - 2 -0.3 




T " 1 "M ; 1 o, n A AO CA/1 

UUr 12-3402-6504 


AQ . Q O .Q /I HO 

03:32:34.02 


T7. /I ^ . CA A 

2 / : 46: 50.4 


7.2 


6.9- 


7 Q. 


— 18.4 


7 Q 7+0.1 
29 - 2 -0.1 


9,13 


T TI~\T771 O. Al O A £H 1 1 

UUr 12-4384-631 1 


AQ . Q Oj . /I Q O A 

03:32:43.84 


—27:46:31.1 


7.2 


6.9- 


7 < 
- / .J 


— 17.9 


OA A _ |-'^ ■ 3 

29 - 9 -0.2 


13 


T TT"\T"?1 O J^f / H -i A 

UUF 12-4256-73 14 


AO O O A O C /" 

03:32:42.56 


/I "7 O 1 /I 

— 27:47:31.4 


7.2 


7.1- 


*7 /I 

- / .4 


— 19.6 


77 7+0.1 

27 - 7 -o.o 


100/)C"7A10 

1,2,3,4,5,7,9,13 


UUF 12-4035-7468 


03:32:40.35 


^*"7 ,1 "7 A f O 

— 27:47:46.8 


7.2 


A *7 

6. /- 


7 £ 

- /.6 


— 17.5 


QA 7+0.4 

30 - 2 -0.3 




UUF 12-3973-6214 


AQ Q "1 OA "71 

03:32:39.73 


— 27:46:21.4 


7.3 


■7 1 
/.l- 


*7 /I 

- / .4 


— 18.4 


on 7+0.1 
29 - 2 -0.1 


1,2,3,4,5,8,9,11 


UDF12-3668-8067 


03:32:36.68 


^—7 1 O A/" "7 

—27:48:06.7 


7.3 


7.0- 


-7.5 


— 18.1 


O A —7 "4" U . O 

29 - 7 -0.2 


9 


UDF12-3708-8()y2 


AO O O O 1 AO 

03:32:37.08 


^ —7 a O AA O 

—27:48:09.2 


7.3 


6.6- 


-7.7 


— 17.5 


on n+0.4 
299 -0.3 




UDF12-4242-6243 


AO O O Af> /I O 

03:32:42.42 


—27:46:24.3 


7.3 


*7 n 
IX)- 


-7.5 


— 18.4 


no n+0.1 

28 - 9 -o.i 


3,5,9 


t rr\r" *1 o o ,1 o 1 "7 i 1 *r 

UDF12-3431-7H5 


AO O O O /I 01 

03:32:34.31 


— * —J j "7 i i — 

—27:47:11.5 


7.3 


7.0- 


-7.5 


— 18.6 


oo c+0.1 




UDF12-3868-5477 


AO OO OO S I > 

03:32:38.68 


—27:45:47.7 


7.3 


6.7- 


-7.8 


— 17.1 


QA Q + 0.4 

30.3+ 03 




UDFl 2-4242-6 1 37 


AO O O .1 O /I O 

03:32:42.42 


—27:46:13.7 


7.3 


6.9- 


-7.7 


— 18.1 


OA O ~ T" (J . 2. 

29 - 2 -0.2 




U Ur 12-41 00-7 216 


AQ 1 AA 

03:32:41.00 


^7.i17.11 £L 

—27:47:21.0 


7.3 


0.5- 


7 P. 
- /.5 


— 17.3 


> 30.6 


13 


U Ur 1 2-4239-6243 


AQ - IT . A 7 Q QA 

03.32.42.339 


— 2/. 46. 24. 3 


7 Q 




7 ^ 
- / .J 


1 O A 

— 1 8.4 


oo 1+0.2 
29 -!-0.2 


Q c n 1 Q 

3,5,9,13 


TTTM-717 /1Q1/1 /:iOC 

UUr 12-4314-6285 


AQ -QT ■ A A It A 

03.32.44.3 14 


— 2/. 46. 28. 5 


7 Q 
/.J 


0.5- 


7 ^ 

- / .J 


1 A 1 

— 19.1 


7Q 1+0.1 
281 -0.1 


1">Q/1C/£"7ni1 1Q 

1,2,3,4,5,6, / , 9,1 1,13 


t rr\c n q q i q /cc/ic 

UUr 12-33 13-6545 


AQ .QO.QQ 1 Q 


— 27:46:54.5 


7.4 


7 

1 .1- 


7 A 
- / .6 


— 18.6 


oo n+0.1 
29 -°-0.1 


1 A C fl 1 1 1Q 

1,4,5,9,1 1,13 


t rr\c n q o o c TC/in 

UUr 12-3885-7540 


AQ .QO.QO CC 

03:32:38.85 


n.yl7.C/l A 

— 27:47:54.0 


7.5 


7 1 

/ . 1 - 


7 7 

— I.I 


— 18.2 


OA Q -j"0 . 2 

29 - 3 -0.2 




TTT~\T7710 QAQ 1 £1 01 

UUr 12-393 1-6181 


AQ . QO.QA Q 1 

03:32:39.31 


— 27:46:18.1 


7.5 


7 1 


7 C 
- /.5 


— 18.2 


on a+0.2 

29.6I .2 


13 


UUr 12-4334-6252 


AQ . Q Oj . /I Q Q /I 

03:32:43.34 


— 27:46:25.2 


7.5 


7 (\ 


-5.U 


— 17.4 


QA A+0-7 

30.6_ . 4 




i jnci o ,-iqao /^O/io 
U L>r 1 2-4308-6242 


AQ -QO-/1 Q AO 

03.32.43.08 


T7. /I 1 

— 2 / .46.24.2 


7 £ 

/.6 


7 1 

/.l- 


7 (1 

- / .9 


— 17.7 


OA £+0.2 

29.6_ n 2 


1 n 1 1 1 q 
1,9,1 1,13 


T ' 1 ' 1 O O O O A "7 A"70 

UDF 12-3880-7072 


AO O "1 OO OA 

03:32:38.80 


n ^ "7 a"7 o 

— 27:47:07.2 


7.7 


7 < 


7 C 

— /.o 


—20.1 


07 7+0.0 

27 - 2 -o.o 


ioo,ir"7Aii 10 

1,2,3,4,5,7,9,11,13 


TTT~\T710 /I O O 1 £CnC 

UUr 1 2-428 1 -6505 


AQ . Q O . /I O O 1 

03:32:42.81 


— 27:46:50.5 


7.7 


7 7 


Q ft 

-5.0 


— 17.9 


OA O ~\~0 ■ 2 

29 - 3 -0.1 




nr\ci o /looo /CQ/lc 

UUr 12-4288-6345 


AQ . Q O . /I O O O 

03:32:42.88 


T7 . /I C -Q /I C 

— 27:46:34.5 


7.7 


7.5- 


-7.9 


— 18.9 


OQ 7+0.1 

28 - 3 -o.i 


1Q/1C£'OA1'1 

1,3,4,5,6,8,9,1 1 


UUF12-4470-6443 


03:32:44.70 


-27:46:44.3 


7.7 


7.6- 


-7.9 


-19.7 


27.8l° ;» 


1,2,3,4,5,9,11,13 


UDF12-4033-8026 


03:32:40.33 


-27:48:02.6 


7.7 


7.5- 


-8.0 


-18.0 


29A±° l 


3,9,13 


UUF12-3722-8061 


03:32:37.22 


-27:48:06.1 


7.7 


7.5- 


-7.9 


-19.2 


28.3+«;i 


1,2,3,4,5,9,11,13 


UDF12-4474-6449 


03:32:44.74 


-27:46:44.9 


7.8 


7.5- 


-7.9 


-18.5 


29.2l° ;| 


1,3,5,13 


UUF12-4240-6550 


03:32:42.40 


-27:46:55.0 


7.8 


7.5- 


-8.1 


-17.9 


29 5+ - 2 


9,13 


UUF12-3939-7040 


03:32:39.39 


-27:47:04.0 


7.8 


7.6- 


-8.1 


-18.4 


28.9t ;} 


13 


UDF12-391 1-6493 


03:32:39.11 


-27:46:49.3 


7.9 


7.5- 


-8.2 


-18.0 


29.61°;* 


9 
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Table Al. Candidate z > 6.5 galaxies in the HUDF continued.. 



Name 


RA(J2000) 


Dec(J2000) 


Zphot 


A 


z 


M1500 


#160 


References 


UDF12-3344-6598 


03:32:33.44 


-27:46:59.8 


7.9 


7.6- 


8.1 


-18.1 


—0.2 


9 


UDF12-3952-7174 


03:32:39.52 


-27:47:17.4 


7.9 


7.8- 


8.0 


-19.3 


28.0+°'i 

— 0.0 


1,2,3,4,5,7,9,11,13 


UDF12-4308-6277 


03:32:43.08 


-27:46:27.7 


8.0 


7.8- 


8.1 


-18.4 


29.3+°-? 

— 0.1 


1,3,5,6,7,8,9,11,13 


UDF12-3780-6001 


03:32:37.80 


-27:46:00.1 


8.1 


7.9- 


8.2 


-18.9 


28.7+!" 

— 0. 1 


1,3,4,5,6,8,9,13 


UDF12-3917-5449 


03:32:39.17 


-27:45:44.9 


8.1 


7.8- 


8.5 


-17.4 


30.1+21 

— 0.3 




UDF12-3762-6011 


03:32:37.62 


-27:46:01.1 


8.1 


7.9- 


8.4 


-17.9 


29.8t°j 




UDF12-3813-5540 


03:32:38.13 


-27:45:54.0 


8.3 


8.2- 


8.5 


-19.1 


28.4±oi 


1,3,4,5,6,8,9,11,13 


UDF12-3763-6015 


03:32:37.63 


-27:46:01.5 


8.3 


8.1- 


8.5 


-18.7 


28.8t°;l 


1,3,4,5,6,9,13 


UDF12-3947-8076 


03:32:39.47 


-27:48:07.6 


8.6 


8.4- 


8.8 


-18.6 


29.0+° \ 


9 


UDF12-3921-6322 


03:32:39.21 


-27:46:32.2 


8.8 


8.6- 


9.2 


-18.0 


29.9±g;l 


13 


UDF12-4344-6547 


03:32:43.44 


-27:46:54.7 


8.8 


8.3- 


9.3 


-17.6 


30.l±g;| 


13 


UDF12-4265-7049 


03:32:42.65 


-27:47:04.9 


9.5 


8.8- 


9.9 


-18.1 


29? +0.2 




UDF12-3954-6285 


03:32:39.54 


-27:46:28.5 


11.9 


11.4- 


12.2 


-19.7 


293 -o 2 i 


12 
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Table A2. Optical and near-infrared photometry for the candidate z > 6.5 galaxies in the HUDF. Detections which are less significant than 2a are listed as 
the appropriate 2a limit. All candidates are detected at less than 2a significance in the -B435, Vqo6 and 1775 filters. 



Name 


^850 


Y 105 


^125 


•A 40 


#160 


UDF12-3999-6197 


> 30.5 


29.81!!;? 


29.7±g;| 


29.9l°j 


29.21°;! 


UDF12-3696-5536 


> 30.5 


99 7+0- 2 
-0.2 


29.8l°-j 


z 9- z to.2 


29.9±g;l 


UDF12-3677-7536 




9S T + 0.0 


28.21";! 


28.31!!;! 


28.21°;! 


UDF12-3897-8116 


> 30.5 


29.8±°J 


?Q 7+0.3 


30 l+°- 4 


29.8l°;3 


UDF12-4120-6561 


> 30.5 


30.0l°j 


30.6l° o j 




> 30.6 


UDF12-3515-7257 


> 30.5 


29.71!!;? 


30.0l°j 


29.9±g;l 


30.1181 


UDF12-3909-6092 


30 4+ ' 7 


29 5+ 1 
z '--'-o.i 


99 7+0-3 
'"'•'-0.2 


30 3+°- 5 


9Q q+0.3 
zy - v -0.2 


UDF12-3865-6041 


30 9+ - 5 
- ,u - -0.4 


29.41J!;! 


29 4+°' 2 
zy/+ -0.2 


99 5+0-2 
zy -0.2 


29.818:2 


UDF12-3702-5534 


> 30.5 


30.0l°| 


30.0l°j 


30.3±°;| 


30.618;! 


UDF12-3922-6148 


> 30.5 


on Q+0.2 


30.21!!;^ 


30.5l ;! 


30 l+°' 4 


UDF12-3736-6245 


> 30.5 


29 5 +o.i 


29.5l° | 


99 5+0.2 

zy "'-0.2 


29 6+°- 2 

zy -°-0.2 


UDF12-4379-6511 


> 30.5 


29.4l!!;i 


Zy - y -0.3 


29.7±g:l 


29.8181 


UDF12-3859-6521 


> 30.5 


29.4l°-J 


29.41^ 


30.0l°j 


29 6+°- 2 

zy -°-0.2 


UDF12-4202-7074 


> 30.5 


29.61!!; 1 , 


29.61- 


90 Q+0-3 


29.5l°; 2 


UDF12-3638-7163 


29.5±°,;| 


28.5i!!;» 


28.51!!;! 


28.41!!;! 


28.61°;! 


UDF12-4254-6481 


> 30.5 


30.0l°;2 


30 3+ ' 5 


30 3+°- 5 


30 6+°- 7 


UDF12-4058-5570 


30 l+°' 5 


29.3l°;l 


29.ll!!;? 


29.4l!!; 2 2 


9q 5+O.2 
zy - -0.2 


UDF12-3858-6150 


> 30.5 


30.3±°;| 


29.9l°| 


30.4l ;« 


30.1181 


UDF12-4186-6322 


30.3t°;« 


29 5+ 1 


29.31°;? 


29.8l°j 


29.418; 2 


UDF12-4144-7041 


> 30.5 


29.9l°| 


90 7+0.3 

Zy -'-0.2 


30.2l°; 5 3 


30.0l°1 


UDF12-4288-6261 


> 30.5 


29.81J!;? 


29.9l°| 


30 l+°- 4 


30.518! 


UDF12-3900-6482 


30.41°;! 


79 1+0-1 
zya -0.1 


29.01J!;! 


28.81!!;! 


28.618;! 


UDF12-4182-6112 


30.51°;! 


29.31J!; 1 , 


29 1+ - 1 


29.01";! 


28.91°;! 


UDF12-4268-7073 


30.11S-J 


29.1±°i 


28.91!!;! 


28.81!!;! 


28.918;! 


UDF12-3734-7192 


> 30.5 


29.31J!;! 


99 0+0.2 


29 5+°- 2 
zy "'-0.2 


9Q 5+0.2 
zy - -0.2 


UDF12-3968-6066 


> 30.5 


30.2l°j 


29.8±g;5 


30.H81 


> 30.6 


UDF12-4219-6278 


29-3t°,i 


28.il!!;!! 


28.ll!!,! 


28.0l!!;° 


28.118;!, 


UDF12-3796-6020 


> 30.5 


29.81!!;? 


9Q 7+0-3 
Zy -'-0.2 


29.8±° l 


9Q Q+0.3 
zy - v -0.2 


UDF12-3675-6447 


30.3l°! 


29.4l!!;i 


29.01!!;! 


29.ll";! 


29.1I8:! 


UDF12-3744-6513 


29.61S-.2 


28 4+°° 
zs -^-o.o 


28.31";! 


28.51!!;! 


28.518;! 


UDF12-4160-7045 


30.01S-J 


28.81!!;! 


28.71";! 


28.71!!;! 


28.818;! 


UDF12-4122-7232 


> 30.5 


29.61!!;! 


29.6l° ; 2 2 


30.0l o ;| 


9Q 7+0-3 
zv -'-0.2 


UDF12-3894-7456 


> 30.5 


90 Q+02 


9Q Q + 0.3 
Zy - y -0.3 


> 30.6 


9Q Q+0.3 
zv - y -0.3 


UDF12-4290-7174 


> 30.5 


30.ll° J 


29.8l°| 


30.2I!!; 4 


30 6+°- 7 

JU -°-0.4 


UDF1 2-4056-6436 


30.2t o j 


28.61!!;! 


28.61!!;! 


28.61!!;! 


28.718;! 


UDF12-4431-6452 


30.0t° o j 


28.7i;;i 


28.61!!;! 


28.61°;! 


28.818! 


UDF12-3958-6565 


zy - v -0.3 


28.41-! 


28.41!!;! 


28.31°;! 


28.418;! 


UDF12-4037-6560 


> 30.5 


29.9±g;l 


29.7l°-J 


30 l+°' 4 


30.1181 


UDF12-4019-6190 


> 30.5 


29.81°;? 


30 1+"' 4 


9Q Q+0.3 
zy - y -0.2 


29.8181 


UDF12-4422-6337 


> 30.5 


30 l+°' 3 


> 30.6 


9q q+0-4 
zy - y -0.3 


30.618;! 
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Table A2. Continued. 



Name 


^850 


5-105 


Jl25 


^140 


#160 


UDF12-4472-6362 


> 30.5 


29.il!!;? 


29.4!8;1 


28.718I 


28.718;! 


UDF12-4263-6416 


> 30.5 


3o.il8:l 


30.1181 


30.41°;° 


30.2l°1 


UDF12-4484-6568 


> 30.5 


29.31°;^ 


29.9!!!;° 


29.4181 


29.6181 


UDF12-3975-7451 


> 30.5 


29.41°,;! 


29.6l° ;l 


29.218;? 


29.4!8;1 


UDF12-3989-6189 


30.51°;* 


29.4l°;i 




29.418:1 


29.5!8;1 


UDF12-3729-6175 


> 30.5 


29.3l°;l 


29.ll!!;! 


29.218? 


28.918;! 


UDF 12-3456-6494 


> 30.5 


29.61"; 1 , 


29.5l° ;l 


29.318I 


29.218;! 


UDF12-4068-6498 


> 30.5 


29.7!!!;? 


29.4l° ;l 


29.4l ;l 


29.6l ;l 


UDF12-3692-6516 


> 30.5 


30.018J 


29.9!8;1 


29.718! 


30.3181 


UDF12-4071-7347 


> 30.5 


29.7!°;! 


30.0181 


29.418;! 


30.1181 


UDF12-4036-8022 


> 30.5 


30.1I8.I 


30.1181 


29.9181 


30.6181 


UDF12-3755-6019 


> 30.5 


29.9!!!;^ 


29.7!8;1 


30.0181 


29-618:1 


UDF12-4256-6566 


29.3l°oJ 


97 o+O.O 


97 i+O.O 

z/ - 1 -o.o 


27.018:8 


26.918:8 


UDF12-4105-7156 


30.41°,;! 


28.7!!!;! 


28.4!!!;! 


28.318;! 


28.418;! 


UDF12-3853-7519 


> 30.5 


29.6!!!;! 


29.4l8;l 


29.418;! 


29.718;! 


UDF12-3825-6566 


> 30.5 


3o.il8:l 


30.2l°1 


29.618;! 


29.818;! 


UDF12-3836-6119 


30.2!°;° 


28.91!!;! 


28.41!!;! 


28.718;! 


28.618;! 


UDF12-3709-6441 


> 30.5 


30.5l°j 


30.5!!!;° 


30.51°;° 


30.2181 


UDF12-3402-6504 


> 30.5 


29.3!!!;! 


28.9!!!;! 


28.918;! 


29.218;! 


UDF12-4384-6311 


> 30.5 


29.9!8;1 


29.518:1 


29.418:1 


29.9181 


UDF12-4256-7314 


30.2!°;° 


28.ii8;° 


27.8!!!;! 


27.818;! 


97 7+0.1 
— 0.0 


UDF12-4035-7468 


> 30.5 


3o.2l8;l 


29.7l8;l 


30.1181 


30.2181 


UDF12-3973-6214 


> 30.5 


29.318;! 


28.81!!;! 


29.318'? 


29.218;! 


UDF12-3668-8067 


> 30.5 


29.6!!!;? 


29.318:1 


29.i!8:? 


9Q 7+0.3 
Zy - ' -0.2 


UDF12-3708-8092 


> 30.5 


3o.2l8;l 


29.9181 


29.818:1 


29.9181 


UDF12-4242-6243 


> 30.5 


29.4!!!;! 


29.018;! 


28.918;! 


28.918;! 


UDF12-3431-7115 


> 30.5 


29.218;! 


28.81!!;! 


28.818;! 


28.5l° : ! 


UDF12-3868-5477 


> 30.5 


30.5181 


30.418J 


30.2l8; 5 3 


30.3181 


UDF12-4242-6137 


> 30.5 


29.718:1 


29.518:1 


29.1181 


29.2!8;1 


UDF12-4100-7216 


> 30.5 


30.518J 


30.0l°1 


30.0l°1 


> 30.6 


UDF12-4239-6243 


30.0l° o j 


29.2!°;! 


29.2!!!;? 


28.918;! 


29.118;! 


UDF12-43 14-6285 


> 30.5 


29.018;! 


28.31!!;! 


28.218'! 


28.118;! 


UDF12-3313-6545 


> 30.5 


29.318;! 


28.91!!;! 


28.618;! 


29.018;! 


UDF12-3885-7540 


> 30.5 


29.8!°;? 


9Q -3+0.2 
-0.2 


29.218'? 


9Q Q+0.2 
-0.2 


UDF12-3931-6181 


> 30.5 


29.8l8;l 


29.ll!!;! 


29.218'? 


29.618;! 


UDF12-4334-6252 


> 30.5 


30.618J 


29.818;! 


29.9l°1 


30.618:1 


UDF12-4308-6242 


> 30.5 


30.3181 


-,-1-0 3 

29.7!8;1 


29.9lg'3 


^r. ^+0 2 
zy -°-0.2 


UDF12-3880-7072 


> 30.5 


28.0l° ° 


27.4!!!;" 


97 -,+0.0 

z/ - -o.o 


97 9+O.O 

z/ - z -o.o 


UDF12-4281-6505 


> 30.5 


3o.5l8:i 


29.6l° ;l 


9Q 7+0.3 
zy -'-0.2 


29.318;? 


UDF12-4288-6345 


> 30.5 


zy - J -o.i 


28.51!!;! 


28.618;! 


28.318'! 


UDF12-4470-6443 


> 30.5 


28.5l°;J 




27.8!8;° 


27.818:8 


UDF12-4033-8026 


> 30.5 


3o.i!8;l 


9Q t+0.2 


29.6!8;1 


29.418:1 


UDF12-3722-8061 


> 30.5 


29.018;! 


28.31!!;! 


28.318'! 


28.318'! 
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Table A2. Continued. 



Name z 850 Y105 ^125 ^140 #160 

2+0.2 2 n +0.l 28 o+O.l 2 n 2 +0.l 
-0.2 zyAJ -0.1 z8 - y -0.1 zy - z -0.1 



UDF 12-4474- 


■6449 


> 


30 


.5 


UDF 12-4240- 


■6550 


> 


30 


.5 


UDF12-3939- 


■7040 


> 


.30 


.5 


UDF12-3911- 


■6493 


> 


.30 


.5 


UDF12-3344- 


■6598 


> 


.30 


.5 


UDF12-3952- 


■7174 


> 


.30 


.5 


UDF12-4308- 


■6277 


> 


.30 


.5 


UDF12-3780- 


■6001 


> 


.30 


.5 


UDF12-3917- 


■5449 


> 


.30 


.5 


UDF12-3762- 


■6011 


> 


.30 


.5 


UDF12-3813- 


■5540 


> 


.30 


.5 


UDF12-3763- 


■6015 


> 


.30 


.5 


UDF12-3947- 


■8076 


> 


.30 


.5 


UDF12-3921- 


■6322 


> 


.30 


.5 


UDF12-4344- 


■6547 


> 


.30 


.5 


UDF12-4265- 


■7049 


> 


.30 


.5 


UDF12-3954- 


■6285 


> 


.30 


.5 



29.%z 2 29.oi;;;j 28.9i;;;i 29 - 2 - 

30.41°;! 29.41°;! 29.71°;! 29.51°;! 

30.11^ 29.1±°;i 29.0l°,;J 28.9±°;1 

30.4t;;; 3 2 29.41";! 29.51°;! 29.61°;! 

on 1+0.3 1O o+0.2 TQ /l +0.2 OQ ^+0.2 



29.3l°;J 28.3Zo:i 28 - 2 lo.i 28 -°lo.o 

3o.2t!!;| 28 - 9 io.i 29 - 4 io.l 29 -3l ;? 

29.81!!;! 28.6±°;1 28.7± ( ( !;^ 28.7±°;1 

> 31.1 29.81!!;! > 30.6 30.ll°;| 

> 31.1 29.51J!;! 29.81!!;! 29.8±g | 

3o.il!!;! 28.61J!; 1 , 28 - 5 -o.i 28 - 4 to:l 

30.6l°j 28.9l°;J 28.9l ( ( !;i 28.81°;} 

31.ll°;I 29.51!!;! 29.0l°;J 29.01°;! 



> 31.1 29.91°;! 29 - 6 -o.2 29 - 9 -o.2 
>31.1 30.0±g;| 30.21°;* 30.1±g;| 

> 31.1 30.6±g:I 29.91°,;! 29 - 7 ta2 

> 31.1 > 30.6 > 30.6 29.3+n ? 
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